arXiv: 1509.00712vl [cond-mat.soft] 2 Sep 2015 


Membrane tension feedback on shape and motility of 

eukaryotic cells 


Benjamin Winkler®', Igor S. Aranson^’*^’^, Falko Ziebert®’^ 

°‘Physikalisches Institut, Albert-Ludwigs-Universitdt, 79104 Freiburg, Germany 
^Materials Science Division, Argonne National Laboratory, 9700 S. Cass Avenue, 

Argonne, IL 60439, USA 

‘^Engineering Sciences and Applied Mathematics, Northwestern University, 2145 
Sheridan Road, Evanston, IL 60202, USA 


Abstract 

In the framework of a phase held model of a single cell crawling on a substrate, 
we investigate how the properties of the cell membrane affect the shape and 
motility of the cell. Since the membrane inhuences the cell dynamics on 
multiple levels and provides a nontrivial feedback, we consider the following 
fundamental interactions: (i) the reduction of the actin polymerization rate 
by membrane tension; (ii) area conservation of the cell’s two-dimensional 
cross-section vs. conservation of its circumference (i.e. membrane inextensi- 
bility); and (iii) the contribution from the membrane’s bending energy to 
the shape and integrity of the cell. As in experiments, we investigate two 
pertinent observables - the cell’s velocity and its aspect ratio. We hnd that 
the most important effect is the feedback of membrane tension on the actin 
polymerization. Bending rigidity has only minor effects, visible mostly in 
dynamic reshaping events, as exemplihed by collisions of the cell with an 
obstacle. 
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1. Introduction 


Motility of cells crawling on substrates attracts substantial interest among 
biologists, physicists, and material scientists alike. Cell motility is a fun¬ 
damental phenomenon that is crucial for a variety of biological processes, 
from morphogenesis to immune response. It is also involved in pathologies 
like cancer growth and metastasis. Like swimming microorganisms, crawling 
motile cells are natural and interesting realizations of active, self-propelled 
systems, displaying self-organized dynamics, flows, as well as intriguing col¬ 
lective effects. Moreover, motile cells and living tissues are inspiring novel 
adaptive materials with intricate properties like active visco-elastic response 
and self-healing. Cellular materials, responding to the topography, elasticity, 
and surface chemistry of the substrate they are in contact with, currently 
inspire microstructured design strategies for cell sorting and guiding. 

The main processes involved in the motion of eukaryotic cells (such as 
keratocytes, hbroblasts or neutrophils) are the following: the generation of 
a propulsive force by actin polymerization against the cell’s membrane, the 
formation of adhesive contact to the substrate to transfer this propulsion 
force and to move forward, and finally, the action of molecular motors in 
determining the cell’s polarity and to retract the rear of the cell (ff. All these 
processes have been modeled in some detail, and models for whole moving 
cells have been recently developed [3, 0, 0, S, 1^. However, there is another 
important player in the game, that has been neglected (or its consequences 
not yet thoroughly studied) in most of the modeling approaches: namely, 
the membrane enclosing the cell. The cell membrane represents a movable 
interface which constitutes an intricate theoretical and numerical problem. 
In addition, membrane tension leads to a global force feedback, affecting the 
propulsion by ratcheting the actin filaments. Moreover, membrane bending 
rigidity may be relevant in some cases, especially for cell collisions with other 
cells or obstacles. 

The first detailed experimental study on the effects of membrane ten¬ 
sion on spreading cells (fibroblasts) dates back no longer than in 2000 [^. 
There, an inverse relation between spreading/lamellipodium extension and 
membrane tension was found: lowering the membrane tension by adding 
detergents (deoxycholic acid) or lipids led to an increased spreading and ex¬ 
tension, while an increase in tension by placing cells in a hypotonic medium 
reduced both effects. The authors concluded that membrane tension may 
constitute a global coupling involved in determining both the cell’s shape 
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and the propulsion dynamics, cf. also the recent reviews [8|, l9[. The effect of 
membrane tension was studied also for neutrophils, both during pseudopod 
formation and for fully developed motion [l^, for spreading fibroblasts (llj |. 
as well as for moving keratocytes [l^. Some of the observed effects include: 
(i) increased membrane tension can cause leukocytes to stop moving [l^; (ii) 
reducing tension can stimulate moving keratocytes to develop several fronts 
[fil ]: ( hi) softening the cell membrane does not affect the velocity of kerato¬ 
cytes [ 12 , ll3|, it only increases the retrograde flow of actin towards the cell’s 
interior. 

Membrane tension has been recently taken into account, for instance, in 
the one-dimensional model for growth cones 1^, as well as the steady state. 
The force balance-based model in [l^ includes also explicit adhesion dy¬ 
namics between the actin cortex and the membrane. Ver y re cently, tension 
gradients and flows inside the membrane were addressed [la, llTJ. However, 


these models do not take shape changes into account, obviously an impor¬ 
tant aspect of the membrane’s feedback. As a result, they can not properly 
describe the onset/cessation of motion. These two important aspects can 
be easily and inherently modeled within t he p hase field approach recently 
developed for motile cells [3, 0, [l^, [l^, 20, 21, 22, 23|, self-propelled active 
droplets (3, 26| and synthetic polymeric capsules j^. Here we include 

and study the most pertinent membrane effects - tension and its feedback on 
polymerization, as well as bending. The study is performed within a simple 
phase held approach for a moving cell. 


2. Phase field model for a crawling cell 

The phase held approach to cell motility has been recently reviewed in 
j^. Instead of modeling the cell’s interface (i.e. the membrane) explicitly, an 
auxiliary held, the phase held p{x,y;t), is introduced. It evaluates to p = 1 
within the cell and to p = 0 outside the cell, with a smooth transition region 
in between describing the ‘smeared’ interface. The simplest implementation 
of the phase held approach is via a scalar order parameter equation 

sp r 

dtp = —, where = J [fip) + Dp{Vpf] dxdy . (1) 

Here /(p) = ^ is a double well potential with minima at p = 0 and 

p = 1 (the two ‘phases’). The phase held free energy Fp in addition includes 
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a surface energy term penalizing interfaces. Equation ([T]) yields 


dtp = DpAp - p{l - p){6 - p) =: X , (2) 


where 5 = | is the ‘pressure difference’ between the two ‘phases’. For 5 = \ 
the free energy of both phases is equal, and hence a planar interface con¬ 
necting states p = 0 and p = 1 is stationary. In case 5 deviates from this 
value, the interface moves either forward or backward, i.e. the cell expands 
or retracts. 

We used this simple framework to model a moving cell by coupling the 
phase held Eq. ([2]) to the polarization held p, describing the averaged local 
orientation of the actin hlaments inside the cell; 


dtp = E>pAp-p(l -p)((5-a|pp -p) -ap- Vp, (3) 

dtp = DpAp - dVp - rf ^p - r 2 ~^(l - p^)p - 7 [(Vp) • p] p . (4) 


In this description, the a-term models the propulsion of the cell’s interface 
by the ratcheting of actin, and the a-term accounts for acto-myosin contrac¬ 
tion. In Eq. dlj), the terms DpAp and —W^p describe dihusion of actin 
and its degradation (depolymerization) in the bulk of the cell, respectively. 
The term —/3Vp describes the creation of actin polarization at the cell mem¬ 
brane (directed normal to the interface) with polymerization rate jS. The 
contribution —r 2 ~^(l — p^)p assures a vanishing polarization outside of the 
cell (where p = 0). Finally, —7 [(Vp) ■ p] p models the front-rear symmetry 
breaking induced by motors. For details we refer to 0,13 


Since motile cells are rather thin (typical lamellipodium thicknesses are 
200 nm) the model is effectively two-dimensional, i.e. height averaged. In 
addition, keratocyte cells are known to preserve their contact area with the 
substrate. To describe this conservation of the cell’s contact area, we intro¬ 
duced the following global constraint 


S = 6v = - + pv [V{t)-Vo] 


( 5 ) 


Here pv is the stiffness of the constraint and the term in brackets is the 
difference between the current area (or 2D volume) V{t) = f p{t) dx dy and 
the prescribed area Vq. Note that, to avoid confusion, in the following area 
always corresponds to the 2D area of the cell’s cross-section (correspond¬ 
ing in a 3D description to the cell’s volume), while the membrane refers to 
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the surface, i.e. circumference, of this cross-section (corresponding in a 3D 
description to the cell’s surface area). 

The position of the interface - which is identified with the cell membrane - 
can be defined in the model by the contour at p = However, this interface is 
not an appropriate description for a cell membrane: it has neither membrane 
tension nor bending energy, but rather an (artificial) wall energy (oc 
that is related to the Ginzburg-Landau-type free energy of the phase field, 
cf. Eq. ([2]). 


Even more important in the context of cell motility is the fact that mem¬ 
brane tension counteracts the polymerization force of the actin hlaments: 
polymerization rate and hence the cell’s velocity decrease as a function of 
the counteracting force, as established theoretically on a single filament level 
by the Brownian ratchet model 29, Although studies of single/few actin 
hlaments polymerizing against a load are very difficult, this effect could also 
been established experimentally (3ll.[3^. 33]. The membrane tension feedback 
on actin polymerization possibly not only leads to a change in the overall ve¬ 
locity of the cell, but also to a global feedback on the actin organization and 
a change in the overall shape of the cell. 


3. Membrane tension as a counteracting force to polymerization 


We will hrst focus on the effect of membrane tension on actin polymer¬ 
ization within the whole cell model described in the last section. To this 
effect, we remove the - artihcial - wall energy of the phase held, and add the 
restoring force of the membrane counteracting polymerization. For simplic¬ 
ity, we keep the simple volume conservation and ignore at hrst the effect of 
tension on the phase held, a limit corresponding to a strongly adhering cell 
that keeps its contact area constant. The effect of tension on the phase held 
is added and studied in the next section. 

The wall energy of the phase held potential can be removed - to leading 
order in the interface width - via addition of the following term to the phase 
held equation 34, dtp = ... + Dpc\Vp\, where c = c(x, y) is the curvature 
of the interface. The latter can be calculated from the local normal unit 
vector, which is determined by the phase held like n{x,y) = -j^, via the 
geometric identity c = —V ■ n. Note that with the given dehnition, the 
normal vector points to the inside of the cell. 

Second, we introduce the membrane tension ({x,y). In the Helfrich pic- 
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ture, the membrane energy reads [35 


/ C\^p\dxdy+- c^\Vp\dxdy. 


( 6 ) 


The hrst term implements the surface area constraint (circumference in 2D), 
where ( is the membrane tension, i.e. the Lagrangian multiplier associated 
with the constraint. The second contribution is the bending energy with 
the curvature c already introduced above and b the corresponding bending 
modulus. The restoring force resulting from this energy has been calculated 
in 361 and reads 


F m.h — 


Ccn — b 


— + t ■ V (t ■ Vc) n + (t ■ VC)t 


|Vp|. 


(7) 


Herein, the first term is the effect of the tension, the second one the contribu¬ 
tion from bending and the last one arises from possible variations in tension 
along the membrane (similar to the Marangoni effect in thermal convection 

0 ). 

Since the polymerization is normal to the membrane, we use only the 
normal contribution F = n ■ F^b- In addition, one can estimate that the 
contribution to the restoring force from bending (oc b) is negligible, see 
, Hence we simply obtain that at the membrane (where |Vp| 
is non-vanishing) the restoring force in normal direction is given hj F = (c: 
the force counteracting polymerization is proportional to the tension and the 
local curvature. 

The simplest way to determine the value of the tension ( is to assume 
that it is related to the overall relative excess circumference of the membrane 
(in 3D; surface area) vial^ 

( = T6A = ^ , where H(f) = (\'\/p\{t) dx dy. ( 8 ) 

^0 J 

Here T is the membrane’s compressibility modulus and Aq is the cell’s cir¬ 
cumference in a reference state. For the latter we chose the stationary, round, 
non-moving state, as it has the smallest circumference. 


Appendix A 


•^Since Vp is nonzero only at the interface, /|Vp|(t) dxdy is (proportional to) the cell’s 
circumference. 
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Finally, we have to account for the feedback the membrane tension pro¬ 
vides on the actin polymerization dynamics. In the simple model, the term 
—/5Vp in the equation for p described that at the position of the membrane 
(where Vp is nonzero), actin is created with constant rate 13 in the normal 
directioqj. Here, we make use of the ratchet-like process of actin polymeriz¬ 
ing against the membrane, where it had been shown that the polymerization 
rate (3 decreases exponentially (in the simplest case) with the force 

13(F) = I3exp (9) 


where ksT is the thermal energy and a is the size of the actin monomer 
2^ 301, which can be absorbed in the constant /q. Note that the parameter 


associated with actin polymerization is not the free polymerization rat^, 
but the rate reduced by the offset tension present in the reference state, 
cf. Eq. dH]). Correspondingly, the propulsion strength, ap, decreases too 
upon an increase in membrane tension, which is due to the smaller amount 
of overall actin polarization p. 

The generalized model then reads 


dtp = DpAp - p{l - p){S - alpl"^ - p) + Dpc\V p\ - ap-Vp, (10) 
dtP = DpAp - Vp - rf ^p - - p^)p - 7 [(Vp) ■ p] p , 

( 11 ) 


with 6 given by Eq. ([5]) and c = —V ■ n = —V ■ fhe curvature as 

described above. For numerical reasons, the normal vector and the curvature 
can only be calculated in a ‘tube’ around the interface (using a threshold value 
for Vp). Upon multiplication with | Vp| the respective terms nevertheless lead 
to smooth contributions at the interface, where they are needed. 

Results for steady moving cells. We investigated the behavior of the 
following quantities that can be easily measured in experiments: the aspect 
ratio of the cell as a measure for the shape change, the cell’s velocity and its 


^Note that Vp is negative when measured from inside the cell, hence the ‘-’sign. 

^ Since /? is not the free polymerization rate, for negative curvature (locally concave 
shape) there could be a slight acceleration of actin polymerization until the free polymer¬ 
ization rate is reached. We neglect this effect here, since it is small and stationary cell 
shapes are almost exclusively convex, by considering only c > 0 in the exponential. In the 
wall energy correction term, however, both curvature signs have to be considered. 
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Figure 1: The aspect ratio’s deviation from round shape, h — 1 (a), the center of mass 
velocity 14 om (b), the relative excess circumference 5A (c), and the membrane tension 
C = T5A (d), for a keratocyte-like cell (black curves; a = 2, /? = 2, cr = 1.2, 7 = 0), 
and for a fibroblast-like cell (red curves; 0 = 3,/? = 1.5, a = 0.9, 7 = 0.3) vs membrane 
compressibility constant T. Remark: the volume is conserved by better then 1% in both 
cases. 

relative excess circumference. Note that the latter is not directly restricted 
via the phase field - only the volume is conserved - but only indirectly via 
the feedback provided on th^olymerization rate. The aspect ratio has been 
characterized as previously |^, by determining the ratio h of the eigenvalues 
of the variance matrix Jjj = f (xi — xf) {xj — x'^)pdx dy, where r'’’ = J rp dx dy 
is the center of mass of the cell. Since a round cell has h = 1, we treat h — 1 
as a measure for the deviation from a circle. 

Figure [U displays a) the aspect ratio’s deviation from the round shape, 
h — 1, b) the center of mass velocity Kom, c) the relative excess circumfer¬ 
ence 6A and d) the membrane tension as a function of the membrane’s 
compressibility modulus T. Shown are results for two different cells, one of 
keratocyte shape (black curves) and one of a more fan-like, hbroblast shape 
(red curves), cf. also figure 2. In case of the first cell type, the aspect ratio 
expectedly decreases by increasing T: at the sides the curvature is highest, 
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Figure 2: The shapes for a) a keratocyte-like cell {a = 2, f3 = 2, a = 1.2, 7 = 0 ) and 
b) a libroblast-like cell {a = ?>, j3 = 1.5, a = 0.9, 7 = 0.3) for small (T = 10) and large 
(T = 200) membrane’s compressibility modulus, cf. also the curves in Fig. [TJ 


and hence the restoring force leads to a more ronnded shape. In contrast, for 
the second cell type, the aspect ratio increases. This is dne to the fan-like, 
triangular shape of the cell, where membrane tension not only reduces the 
extension normal to the direction of motion, but also in direction of mo¬ 
tion. The velocity is not substantially affected (in the 10 — 20% range), in 


accordance with experiments [l^, [12 


The relative excess surface decreases in both cases, cf. Fig. [TJ:;), which 
shows that in the limit of high T, the surface (circumference in 2D) should 
become rather well conserved, even without explicit inclusion into the phase 
held equation. Although the relative excess surface decreases, the overall 
tension C, = T5A increases, cf. Fig. [TJi), as it should. Experimentally, also 
the membrane tension is accessible by pulling membrane tethers from moving 
cells, as has been recently studied in |l2| . 

Figure [2] shows the respective shape changes, comparing low and high 
membrane extensibility moduli (end hence tension values) for the two kinds 
of cells. 

Arrest of cell motion by increased tension. Next we have stud- 
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Figure 3: Evolution of the aspect ratio’s deviation from round shape, /i — 1, and center 
of mass velocity 14om for mildly driven cells (parameters: a = 1.3, /? = 1.3, cr = 1.4, 
7 = 0). The tension modulus T was increased step by step and the velocity and aspect 
ratio measured after the cell had equilibrated to its steady state. At a certain critical 
value Tc the cell stops. The snapshots at the left and right show the cell’s shapes for small 
(negligible) tension and for the highest tension value where the cell still is able to move. 


ied whether the onset/cessation of motion is affected by membrane tension. 
Mildly driven cells (i.e. cells with not too large values of the propulsion pa¬ 
rameter a, as well as a and 7 ) can indeed be stopped by increasing the 
compresibility modulus T and hence the tension. This effect has been seen 
in several experiments, e.g. in [l^ for leukocytes, and is captured by our 
model as shown in Fig. [3l Interestingly, while without membrane tension 
both velocity and aspect ratio exhibit a jump at the arrest of motion, cf. j^, 
with tension present the aspect ratio decreases continuously and the cell be¬ 
comes even slightly stretched in the direction of motion, before hnally being 
stopped for even higher tension. The reason is that the propulsion is pre¬ 
dominantly due to the front-rear asymmetry. Increasing T strongly affects 
the sides (i.e. the aspect ratio decreases), but the feedback on polymerization 
is small at the front (and typically also at the rear) since curvature is small 
there. Consequently, the jump in the velocity is decoupled from the jump in 
the aspect ratio, due to tensioij^. 

Experiments indicated a direct relation between membrane tension and 


® Note that the study in Ref. @ indicated that there is not a simple relation between 
the aspect ratio and the center of mass velocity, e.g. when varying the cell’s contact area 
(volume in 3D); here the same is true even for fixed contact area, but for different values 
of tension. Hence again, the aspect ratio is not necessarily related in a simple way to the 
cell’s velocity. 
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the actin pushing force: reducing the latter results in a decrease in membrane 
tension, as found by treating keratocytes with blebbistatin and jasplakinoli- 
dine, leading to a rapid cessation of the actin assembly jl^. Ref. [l^ found 
the same, namely that membrane tension increases upon stronger leading- 
edge protrusion. This relation can be directly inferred from our model: 
namely, it is the propulsion by actin oc a (and caused by actin polymer¬ 
ization oc /3) that induces an increase in the excess area 5A (even in the 
non-motile case) and hence in tension C, = T 6 A/A 0 , cf. Eq. ([ 8 ]). 

Another interesting experiment performed in Ref. [l^ was the fusion of 
a moving keratocyte with a giant unilammelar vesicle (GUV) to increase 
the available membrane area. The cross-sections of the cells were found to 
become larger (i.e. volume entered from the third dimension) with higher as¬ 
pect ratios, while the velocity and tension remained the same. This behavior 
too is in accordance with our model: after fusing a cell with a GUV, both 
the contact area Vq and the circumference Aq will increase due to spreading. 
We already studied the behavior of velocity and aspect ratio as a function of 
the cell’s contact area in j^, and found that (keeping all other parameters 
hxed) for not too large cells the aspect ratio increases with area while the 
velocity remained practically unchanged. Membrane tension will not change 
this scenario due to the concomitant increase in Aq (the restoring force may 
decrease a little, since it is proportional to curvature, but the leading edge is 
rather strai ght anyways). 

In Ref. [l2| it was also observed that an increase in cellular adhesion 
to the underlying substrate increased the membrane tension. Again, this is 
consistent with our model since the propulsion force is proportional to the 
number of adhesive ligands A, a ~ aoA, see Ref. [l^ for a generalization 
of the model including explicit adhesion dynamics. Gonsequently, upon in¬ 
creased adhesion the cell can spread more efficiently, thereby increasing the 
relative excess area and consequently the tension 5AT. 


Finally, Ref. [12| found that a decrease of myosin contraction leads to 


higher tension. This is the only trend not (yet) captured in our model: if the 
parameters associated with the activity of motors {a or 7 ) are decreased, we 
do obtain more round shapes, but this is not due to an increase in tension. 
The reason for this discrepancy might well be that the implementation of the 
action of motors is still too oversimplified: we neither implemented explicit 
motors, nor tensorial active stresses, nor the contractile bundle at the rear 
present in keratocytes. 
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4. Contact area vs. contour conservation 

There arises the question, which quantity should be conserved in an ef¬ 
fective 2D model of a cell. For a 3D cell this is rather clear: the volume is 
conserved since the cytoplasm is incompressible, and the membrane area too 
since the membrane’s compressibility modulus is very high, implying almost 
perfect inextensibility (note, however that there are cells with membrane folds 
to buffer surface area, see also below). In contrast, for a height-averaged 2D 
model as ours, both the contact area could vary - e.g. the cell retracts to the 
third dimension by reducing its spreading - and the circumference. 

One can describe this effect on a phenomenological level by considering 
the phase field parameter 6 to depend on both the contact area and the 
length of the circumference: 

S = - + fiy [V{t) ~ bo] + f^AT [A{t) — Aq] . (12) 

Here, we tuned /iy and in such a way that for a moderate value of the 
compressibility modulus, T = 50, both contributions are of same order (for 
a cell of specific size; we typically used cell’s of radius tq = 15). Hence T = 0 
corresponds to pure contact area conservation, while the limit T S> 50 leads 
to a dominating conservation of the circumference. 

We investigated the model described by Eqs. ffTOj) and flTT]) . with the pa¬ 
rameter 6 replaced by Eq. ([12]). Figure ID displays the aspect ratio, the center 
of mass velocity, the contact area dV, and the circumference 6 A, as a function 
of the membrane’s compressibility modulus T. Expectantly, with combined 
area-circumference conservation, the circumference is better conserved at the 
cost of the contact area conservation [see panels c) and d)]. Figures |4^) and 
b) show that the effect of tension on both the aspect ratio and the velocity are 
much smaller than when only area conservation constraint is imposed. Hence 
interestingly, although one would naively expect a decrease of speed from the 
exponential decrease of polymerization with tension, the speed remains al¬ 
most unchanged due to the global shape dynamics: one can not conclude the 
overall speed just from the local polymerization rate. In case of only con¬ 
tact area conservation, the velocity decreased by 10-20%, cf. also Fig. [T] In 
contrast, if the circumference conservation is important, the overall decrease 
is below 5% and probably not even measurable in experiments. Note that 
keratocytes are known to have only transient, weak adhesion and high mem¬ 
brane tension [l^, hence circumference conservation should be important. 
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Figure 4: Shown are the aspect ratio’s deviation from round shape h — 1 (a), the center 
of mass velocity Kom (b), the excess contact area 6V (c), and the relative excess circum¬ 
ference 5A (d) for two cell types as a function of the membrane’s compressibility modulus 
T: for a keratocyte-like cell (a = 2, /3 = 2, cr = 1.2, 7 = 0) with area conservation (black 
solid) and with combined area-circumference conservation (black dashed), as well as for a 
fibroblast-like cell (a = 3, /3 = 1.5, a = 0.9, 7 = 0.3), again with area conservation (red 
solid) and with combined area-circumference conservation (red dashed). Note that the 
solid curves are identical to those in Fig. [H 


Indeed both experiments 13|, ll2(] did not see an effect of membrane tension 


on the cell’s velocity. We are not aware of such experiments on hbroblasts, 
but given their higher adhesion and lower tension, the study in the previous 
section suggests that their velocity might be affected in the 10-20% range. 

Thus, our study indicates that the exchange of contact area and mem¬ 
brane circumference with the third dimension is crucial for understanding the 
influence of membrane tension on the shape and speed of cells. Obviously, 
this problem can only be properly addressed in truly three-dimensional mod¬ 
els which are underway 0 - Note that the membrane’s curvature in this 


third direction is typically much higher than the in-plane one (due to the 


thinness of the lamellipodium) and might be the dominating curvature 38 


Another complication arises due to the existence of plasma membrane reser¬ 
voirs. While these reservoirs seem not to be present in keratocytes 0 , they 
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are known to exist for neutrophils (where they are important for phagocyto¬ 


sis 


39| . i.e. the uptake of micron sized objects), as well as in hbroblasts. The 


percentage of area buffered/stored can be between 1% for strongly adhering, 
up to 10% for spreading hbroblasts [ll|, 4^ . 


5. Effect of bending rigidity &; perturbing moving cells 

The bending rigidity contribution to the force opposing actin polymer¬ 
ization is typically negligible compared to the contribution due to tension 
(cf. Section [3] and Appendix). Nevertheless, bending may affect the shape of 
the cell directly when its contribution to the phase held equation becomes 
relevant. From the typical 2D radii of curvature (of order 20 pm) and typical 
membrane rigidities, one expects only moderate ehects for stationary mov¬ 
ing cells - mostly at the wings where the curvature is the highest. However, 
bending may well be relevant for the dynamics of cells and for the response 
to external perturbations, e.g. via modulations of the substrate properties as 


studied in [19|, 12^ or by obstacles 


For vesicles, the bending energy can be introduced similar as the phase 
held energy in Eq. ([T]), by considering the already introduced function A = 
DnAp— f'(p). Helfrich-Willmore theory then implies the bending energy 
to be lial 


E,= 


6^2 ^ 
Ep\/T^p 


}? dx dy. 


(13) 


In the sharp interface limit, see Ref. (42|, this expression tends tob f Cf^dx dy, 
where Cm is the mean curvature and b is the bending modulus. However, 
this approach can only be used in the ‘advected held’ case, i.e. if the shape 
of the phase held p across the interface stays (close to) a tanh-prohle all the 
time. This is the case for vesicles in Stokes how, and also in the active gel hy¬ 


drodynamics approach recently studied 25(] , since there are no forces normal 
to the interface. In the model for a moving cell considered here, however, 
the active terms (proportional to a and a) are normal forces and deform 
the interface away from the radial tanh-shape. We checked the overall bend¬ 
ing energy using Eq. (|T3l) for stationary round cells. We obtained that for 
a > 0.5 - already well below values needed to induce motility - the bending 


energy not only increased in value, but, in addition, did not display anymore 
the correct behavior E;, oc 4. Hence in Eq. (irsj) one would have to correct 
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for the terms perturbing the radial profile of the phase held, which is neither 
a straightforward task nor very intuitive. 


We therefore applied a purely 2D formulation for the bending energy [36 
given by Eq. ([6]), i.e. used 


Eb = - J c^|Vp| dxdy, 

leading to the following contribution in the phase helc0 


d,p=... + b<>- + t-V{t-Vc)}\Vp\. 


(14) 


(15) 


As expected, for steady moving cells the effect of bending is not very 
noticeable, except for very low tension values and rather strongly elongated 
cells. Figures [5^) and b) show steady moving cells with T = 10, and b = 0 
and 1, respectively. One can see that bending leads to a rounding of the 
wings and slightly also of the back. 

We also investigated collisions of cells with rigid round obstacles, to test 
whether perturbations externally imposing a curvature lead to stronger ef¬ 
fects. For this purpose, the obstacle was implemented by a second stationary 
phase held po{x,y) (that had been relaxed towards a radial tanh-prohle of 
equal width as the stationary cell) and coupled to the cell via a steric inter¬ 
action energy (as developed previously for multiple cells (23[) 


A'steric Pi 


kplp. 


( 16 ) 


For the interaction strength we chose k = 10, typically, to prevent overlap. 
Figures [5hi)-d) shows a collision with an obstacle of radius Ro = lO/rm for 
a cell with low tension, T = 10, and high bending rigidity, 6 = 1. For 
comparison. Figures [5^)-h) display a cell with the same bending rigidity but 
higher tension T = 50. One can clearly see that the cell with low tension 
becomes very deformed in spite of its bending rigidity, while the second one 
much less. The second cell relaxes rapidly to its initial shape, albeit dehected 
in its direction of motion, cf. also Suppl. Movies 1 & 2. 


^ Since the bending contributions scale with the radius of curvature like 1/i?^ and 
rapidly grow towards the inside of the cell, for numerical reasons one has to decrease Dp 
slightly (increasing the sharpness of the interface) and also the tube width around the 
interface where the terms are calculated. 


15 








Figure 5: Sequences of snapshots illustrating the collision between a moving cell and an 
obstacle. Panel a) shows a steady moving cell with low tension (T = 10) and no bending. 
Panels b)-d) display the collision of a cell with the same tension and bending rigidity 6=1 
with a round obstacle. From the shape of the steady moving cell b) in comparison to a) 
one can infer that bending leads to a slight rounding of the wings and the back. For such 
a ‘floppy’ cell of low tension, the collision leads to a strong deformation of the wing hitting 
the obstacle. Parts e)-h) show a collision of a cell with higher tension T = 50 and same 
bending rigidity, with an obstacle. The cell is much less deformed, recovers its initial shape 
more rapidly and is slightly deflected from its initial direction of motion. The collision 
events have been traced in the moving frame of the cell. Other parameters: a = P = 1.3, 
a = 0.6, 7 = 0.4; Dp = 0.5 in place of Up = 1 which was used in the other figures. 


Finally, we also studied cells encoutering a narrow channel, a situation 
increasingly studied experimentally 0, cf. also H- Expectantly, for large 
enough channel widths the cells pass, while narrow channels stop them at the 
entrance. Suppl. Movies 3 & 4 focus on the intermediate behavior: Movie 
3 shows a cell with low tension, T = 10, that is able to go through the 
channel albeit being slowed down. Movie 4, in contrast, shows a cell with 
high tension, T = 50, encountering the same channel. It is less deformable 
and hence does not pass. 


6. Conclusions 

We investigated the consequences of a variety of physical mechanisms 
associated with the cell membrane, such as the feedback of membrane ten¬ 
sion on actin polymerization, contact area vs. circumference conservation. 
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and membrane bending stiffness, on the shape and motility of moving cells. 
These questions are important in the context of cell motility, since especially 
the membrane tension acts as a global mechanical feedback that may con¬ 
stitute a conduit for rapid - note that tension relaxes on the time scale of 
few milliseconds - information transfer across the cell. In addition, phys¬ 
ical perturbations of cells, by e.g. collisions or substrate modulations, are 
increasingly investigated. The overall conclusion is that the membrane plays 
an important role in preserving the integrity of the cell as well as for its 
dynamical response. 

We have found that the dominant effect of the membrane on the cell’s 
motility is its feedback via tension on the ratcheting of actin. This is in accor¬ 
dance with the hnding by Lieber et al. [l^ that for keratocytes, membrane 
tension is dominated by the cytoskeletal forces. Since keratocytes have rather 
high tension compared to other moving cells, the most appropriate model for 
this cell type is the one with combined contact area and circumference con¬ 
servation, developed in section 01 Accordingly, bending only leads to small 
corrections, mostly during cellular reshaping events. 

Interesting future aspects of membrane effects include generalizations to 
3D models of cells. There, the high out-of-plane curvature at the tip of the 


lamellipodium may dominate the behavior [38j . Another important aspect is 


related to the fact that the tension originates not only from the lipid bilayer 
membrane, but also from the actin cortex directly underneath, adhering to 
the membrane. A loss of the membrane-cytoskeleton adhesion then can lead 
to blebbing jd^, constituting another means of motility, which has not yet 
been modeled in a whole cell, dynamical model. 
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Appendix A. Parameters and estimate of tension vs. bending 

The typical scales in our model are chosen to be seconds for time, microns 
for length and piconewtons for force (see [sj for details). 
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The force given in Eq. ([7]) was derived for vesicles and is a volume 
force (entering the Stokes equation), i.e. it has units of N/m^. To obtain 
the force acting on an actin hlament, we hence have to multiply by the 
volume of an actin monomer, roughly = (10 nm)^. In the feedback on the 


polymerization rate, in the expression exp one also has a = 10 nm and 

ksT ~ 4pNnm. 

The membrane’s compressibility modulus is of order T = O.lJ/m^ 
leading to 10® pN//im. For numerical reasons we varied T = 0 — 500. The 
tension C should be of the same order j^; it has been recently measured in 
[l^ by pulling membrane tethers from moving keratocytes, and was found 
to be of the order of 250 pN//im. The force due to membrane tension can 
be hence estimated to be (Cc|Vp|a® ~ 0.01 — 1 pN, depending on whether 
one uses the typical phase held interface width (~ 1 /rm) or rather a real 
membrane thickness (5nm). Consequently, the prefactor /o in exp(—/ qCc) 
occurring in Eq. (ITT]) should be of order 0.01 — 1, we typically used 1. 

In contrast, for the bending rigidity contribution in Eq. ([7|) we obtain 
5c®|Vp|a® = 10“® — 10“^pN, using a typical membrane bending modulus of 
order oi b = 100 ksT = 0.4pNyum. In 1^ also the bending rigidity has been 
directly measured for keratocytes and found to be slightly lower, 0.14pN/im. 
Thus, for the membrane’s feedback on the polymerization rate, bending is 
completely negligible vs. tension. In the phase held equation, however, the 
bending contribution enters directly i.e. is of order 0.14 — 0.4 in our units, 
we typically used 6 = 0 — 1. 
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